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Abstract. A system of globally coupled rotors is studied in a unified framework 
of microcanonical and canonical ensembles. We consider the Fokker-Planck equation 
governing the time evolution of the system, and examine various stationary as well as 
non-stationary solutions. The canonical distribution, describing equilibrium, provides 
a stationary solution also in the microcanonical ensemble, which leads to order in a 
system with ferromagnetic coupling at low temperatures. On the other hand, the 
microcanonical ensemble admits additional stationary and non-stationary solutions; 
the latter allows dynamical order, characterized by multiple degrees of clustering, for 
both ferromagnetic and antiferromagnetic interactions. We present a detailed stability 
analysis of these solutions: In a ferromagnetic system, the canonical distribution 
is observed stable down to a certain temperature, which tends to get lower as the 
number of Fourier components of the perturbed distribution is increased in the analysis. 
The non-stationary solution remains neutrally stable below the critical temperature, 
indicating inequivalence between the two ensembles. For antiferromagnetic systems, 
all solutions are found to be neutrally stable at all temperatures, suggesting that 
dynamical ordering is relatively easy to observe at low temperatures compared with 
ferromagnetic systems. 
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1. Introduction 

The system of sinusoidally coupled oscillators serves as a prototype model describing 
various oscillatory phenomena in nature. When the coupling is short-ranged, i.e., 
between nearest neighbors, the oscillator system describes an array of Josephson 
junctions, which has been a subject of extensive studies On the other hand, there 
are also many systems with long-range couplings in physics and biology. Physiological 
rhythmic processes may be examples of the latter, which may be modelled as a system of 
coupled oscillators with the range of coupling being varied, where phase synchronization 
of the system is an important issue to be understood Pj. Physical examples are 
diverse, ranging from self-gravitating and plasma systems, where the long-range nature 
of the gravitational or Coulombic interaction gives rise to difficulty in understanding 
the systems. A system of globally coupled rotors has thus been proposed and studied to 
simulate those systems [3]. Here the interaction range is infinite, with the strength scaled 
with the system size, making the system of the mean-field character and amenable to 
analytical treatment. In spite of the mean-field nature, however, the system has turned 
out to exhibit rich features in dynamical and statistical properties. 

In the canonical ensemble one can find an analytic solution and the system with 
the ferromagnetic interaction undergoes an equilibrium phase transition at a finite 
critical temperature, whereas there is no phase transition for the antiferromagnetic 
interaction. On the other hand, direct simulations in the microcanonical ensemble 
reveal some interesting features with remarkable differences with the nature of the 
interaction. Specifically, for the ferromagnetic interaction, the system displays extremely 
slow relaxation towards the thermodynamic equilibrium. This slow relaxation, dubbed 
quasi-stationarity, does not coincide with predictions in the canonical ensemble, and 
thus suggestion has been made that there may exist inequivalence between canonical 
and microcanonical ensembles. Such quasi-stationarity is observed to survive well 
below the equilibrium critical temperature and hence has attracted much attention 
[3 m 13 El 13 El El Uni m] , together with some controversy [Zj . In the regime showing 
quasi-stationarity it has also been reported that the system exhibits aging effects and 
glassy behavior [71 |H]. For the antiferromagnetic interaction the system exhibits a 
different type of coherent motion at low temperatures, again only in the microcanonical 
ensemble The rotors move in two groups, called the bi-cluster, for a long time, 

which is explained in terms of the statistical equilibrium of the effective Hamiltonian 
obtained after averaging out fast variables. 

In a recent work we have employed a novel approach that treats the system in a 
unified framework of microcanonical and canonical ensembles fSj- Starting from the set 
of Langevin equations describing dissipative dynamics of a system (canonical ensemble) 
and the corresponding Fokker-Planck equation (FPE), we have pointed out that the 
nondissipative Hamiltonian dynamics (microcanonical ensemble) may be described as a 
limiting case of the vanishing damping coefficient. Thereupon we have been able to find 
a class of solution for the incoherent phase depending on the ensemble, some of which are 
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neutrally stable even below the equilibrium critical temperature. This neutral stability 
has then been suggested to be a plausible physical explanation as to the origin of the 
quasi- st at ionarity observed in numerical experiments. In this paper we further extend 
the stability analysis of the previous work to the ferromagnetic coherent phase (with 
thermodynamic order) and to the system with the antiferromagnetic interaction. For 
the latter, we attempt to provide an alternative view of the bi-cluster phase observed in 
the antiferromagnetic system, as dynamical order allowed by the rotating solution of the 
FPE. This rotating solution is found to be neutrally stable down to zero temperature. 
Furthermore, the rotating solution can give rise to any degree of clustering, if the initial 
condition is appropriately chosen, in addition to bi-clustering. It would thus be of 
interest to probe such multi-cluster motions as tri-clustering, etc., by means of numerical 
simulations. 

This paper is organized as follows: In Sec. El we describe how the system of globally 
coupled rotors can be treated in a unified framework from the set of Langevin equations 
and the corresponding FPE. Various solutions of the FPE are given in Sec. El It is 
shown that multi-cluster solutions emerge, manifesting dynamical order for the non- 
stationary rotating solution of the FPE. Section (His devoted to the stability analysis of 
the stationary solutions, with emphasis on the ferromagnetically coherent phase (with 
single cluster motion or thermodynamic order). The stability analysis of the non- 
stationary solution is presented in Sec. El with a special focus on the antiferromagnetic 
case. Finally, a brief summary is given in Sec. 



2. System of Coupled Rotors 

We consider a system of classical rotors, each of which is described by its phase 
angle and coupled sinusoidally to others. The dynamics of the coupled rotor system is 
governed by the set of equations of motion for the phase (pi {i = 1,...,N) of the ith 
rotor: 

Mcbi + J2 J^j M(t>i - = 0, (1) 

3 

where M is the inertia of each rotor and Jij represents the coupling strength between 
rotors i and j. With the introduction of the canonical momentum pi = M(j)i, the above 
equations are transformed into a set of canonical equations: 

with the A^-particle Hamiltonian 

'^N = ^l^-^JijCOs{(t)i- (pj), (3) 
i i<j 

on which the microcanonical description is based. 
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On the other hand, in the canonical description the system is in contact with a heat 
reservoir of temperature T and described, in a most general way, by the set of Langevin 
equations: 

M^i + r(j)i + J2 Jij sin(</>^ -<l)j) = r]i, (4) 
j 

where F is the damping coefficient and the Gaussian white noise rii{t) is characterized 
by the average (rjiit)) = and the correlation {Vi{t)rjj{t')) = 2TTSijS(t—t'). To derive 
the corresponding FPE, we write the equations of motion in the form 

r 

It is then straightforward to derive, via the standard procedure JHI, the FPE for the 
probability distribution P{(f)i,pi,t): 

dP _ ^ dP ^ d 

P (6) 



X 



F d 
-rPi + X] Jij sin(0i - (t>j) + FT- 



One may also derive the FPE for the Hamiltonian dynamics, which just reads Eq. (jH)) 
with F = 0. While reflecting that Eq. (jH) with F set equal to zero reduces to Eq. ((T)), 
this suggests that Eq. © should provide the starting point for both descriptions: the 
microcanonical one (F = 0) and the canonical one (F 7^ 0). In particular, the stationary 
solution of Eq. is given by the canonical distribution P^^\(j)i^pi) oc e^^'^/"^, 
describing equilibrium, with the very Hamiltonian in Eq. Q regardless of F being 
zero or not. Note, however, that unlike the canonical ensemble where T represents 
the given temperature, in the microcanonical ensemble T still remains as an arbitrary 
parameter. In the latter, one may adjust T to the average kinetic energy, which 
allows the interpretation of T as the temperature. This prescription thus establishes 
correspondence between the two ensembles. Note also that in the zero-temperature 
limit (T 0), Eq. reduces to the Caldirola-Kanai Hamiltonian dynamics jl6j . 
which needs external driving to have a nontrivial stationary state. 

In order to measure a variety of coherence in the system, we conveniently introduce 
the generalized order parameter A*^^) defined by 

1 ^ 

-^e'^<^^ = A^e'^^ (7) 

i 

Apart from the global phase non- vanishing values of the order parameter A^^-* 
imply that rotors move as clusters, since rotors separated with phase angle 27r/£ make 
contributions to A^^^. It thus can be used as a measure of the distribution of rotors, 
particularly, the degree of clustering. For instance, a non-vanishing value for 1 = 1 
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corresponds to the emergence of a mono-cluster (or magnetization), that for i = 2 
corresponds to bi-cluster formation (with separation of vr), and so on. Note that the 
i = 2 case may be regarded as the analogue of staggered magnetization in the short- 
ranged model. 



3. Stationary and Non-Stationary Solutions 

In the infinite-range limit (J^- = J/N with oo), we use Eq. (j?)) for i = 1, 

A« = ^5^e^('^-^^), (8) 

i 

and decouple the set of the equations of motion into a single-particle equation 

M0; + + JA(i) sin(0i - 9,) = rii, (9) 

satisfied by all rotors. Henceforth we therefore drop the rotor index i in Eq. 0, which 
leads to the standard FPE for the single-rotor probability distribution P{(f),p,t): 

dP p dP , wl^ . ^,dP 
— = - —— + JA(^) sm((h - 9i) — 
dt M (90 ^ dp 

P (10) 



dp 



M dp 



In the absence of damping (F = 0), this reduces to the FPE for the micro canonical 
ensemble: 

dP p dP . ,^ ^ ,dP 

which is also referred to as the Vlasov equation in some literature [SlElEl- In terms 
of this probability distribution, the generalized order parameter is defined to be 

3.1. Stationary Solutions 

For the sake of completeness, we briefiy review the results for stationary solutions of 
the FPE [ini. As pointed out for the general case, both Eqs. (fTUj) and (fTTj) support the 
same stationary {dP/dt = 0) solution: 

P(°H0,p) = ^e-«/^ (13) 
with the single-particle Hamiltonian 

2 

n = ^-JA^'^cos{<P-9^), (14) 

where the overall phase 9i manifests the global U{1) symmetry. It is thus expected that 
both ensembles exhibit the same equilibrium behavior. 

One, however, should recall again that here T is given in Eq. ()10j) (for the canonical 
ensemble) but remains arbitrary in Eq. (jlip (for the microcanonical ensemble). In the 
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microcanonical ensemble the temperature should be defined as a measure of the average 
kinetic energy according to (jp) /2M = T/2. The partition function is determined by 
normalization: 

'dpj^e-^/^. (15) 

For later use, we first describe some equilibrium properties of the globally coupled 
rotors [3] through the use of the single-particle model. Defining x = J/SS^^/T and 
making use of the expansion 



gXCOs(0-ei) 



^ /„(x)e^"(<^-^^) (16) 



with In{x) being the modified Bessel function of the n-th order, we evaluate the partition 
function as 



Z = V27iMTIo{x). (17) 

We emphasize again that this approach based on the FPE provides a unified 
description of microcanonical and canonical ensembles and both ensembles generate the 
same equilibrium behavior, determined by the same distribution P^^\(j)^p). Namely, in 
both ensembles the generalized order parameter in equilibrium is given by 



With the expansion in Eq. ()lfi|l and integration over 0, the order parameter reads jU! 

Note here that A*^^) has an explicit dependence on the coherence order parameter A*^^) 
through X = JA^^^/T. For i = 1, describing the emergence of coherence (the mono- 
cluster as thermodynamic order), Eq. (|THjl becomes an equation to be solved self- 
consistently: 

J Io{x) 

This self-consistency equation determines whether the system exhibits coherence: The 
ordered phase (A*^^^ ^ 0) emerges when T/J is smaller than the slope of Ii{x)/Io{x) 
at X = 0, which is 1/2. Accordingly, the ferromagnetic system (J > 0) undergoes a 
phase transition at the critical temperature = J/2. In the case of antiferromagnetic 
coupling (J < 0), on the other hand, Eq. ^lU^ becomes —Tx/\J\ = Ii{x)/Iq{x), leading 
to the only solution x = 0. It is thus concluded that the antiferromagnetic system has 
no phase transition at finite temperatures (no mono-cluster). It is obvious in Eq. (fTn|) 
that A^^^ for higher values of i can assume nonzero values only for A*^^) ^ 0; this implies 
that only the ferromagnetic system can develop all degrees of clustering below Tc. This 
is not surprising since the mono-cluster phase has a 2it symmetry and therefore invariant 
under any rotations of multiples of 27r, which in turn gives rise to nonzero A^^\ 
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For A(^) = 0, describing the incoherent phase, the single-particle Hamiltonian (jl4j) 
has only the kinetic energy term, thus reducing the canonical distribution P^°^((/),p) 
to the Maxwell distribution for both ensembles. Unlike Eq. (fTIUl . however, Eq. 
flll|) . the FPE in the microcanonical ensemble, allows an extra solution of the form 
P''°-*(0,p) = /o(p), an arbitrary function of p without (/)-dependence, including the 
Maxwell distribution as a special case J^]- The only constraint is the normalization, 
and the distribution P(°^(0,p) uniform in guarantees A*-^) = 0. As a result, A^^^ 
vanishes for all values of i as well, and no multi-clustering is allowed by this type of 
stationary solution present only in the microcanonical ensemble. 



3.2. Rotating Solutions 

In addition to the stationary solutions presented above, the FPE in the microcanonical 
ensemble also carries non-stationary solutions which have some significance for the 
antiferromagnetic system. For A*^^-* = 0, Eq. (jllj) becomes 

dt M del)' ^ ' 

which has a solution of the general form P^^\(l),p, t) = u{(j)—pt/M,p). This is a rotating 
solution in the sense that the phase grows continuously with time with a continuous 
frequency spectrum {u oc p/M). Requiring periodicity in 0, we write 

P(°)(0,p,t) = 5^e^'=(<^-^*)Pfc(p), (22) 

k 

where Fk{p) is an arbitrary function of p satisfying P±i(p) = due to the condition 
A*^^) = 0. The generalized order parameter for this solutions is computed according to 



A^e^^^ = j dp j g^e^^V'^^^-i^^^P.b) 



k 

dpe+^'^'F.eip), (23) 

which shows that the higher-order moment A*^^) in general does not vanish unless 
F-ilp) = 0. Thus far there is no difference between the ferromagnetic and the 
antiferromagnetic couplings. As will be shown later, however, the stability of the 
rotating solution differs substantially, depending on the nature of the interaction. The 
rotating solution exists only for A*^^) = 0, regardless of whether the system is in 
equilibrium or not. While such an incoherent phase appears only at high temperatures 
in the ferromagnetic case, remains always zero at all temperature ranges in the 
antiferromagnetic case. Moreover, the rotation frequency gets higher as the order of the 
moment increases. This suggests that at low temperatures where thermal fluctuations 
are small, the phases with non-vanishing higher moments (high degrees of clustering) 
are easier to observe in the antiferromagnetic system than in the ferromagnetic one. In 
fact, this is precisely what has been seen in recent numerical simulations, which reported 
the bi-cluster phase in the antiferromagnetic system at very low temperatures 
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The bi-cluster state, with two clusters separated by angle vr, may be obtained with 
suitable choices of Fk{p) in Eq. ()22p. For example, with the choice F2k{p) = F{p) and 
F2k+i{p) = 0, we obtain [T3j 



P'^'^\4>,P,t) = ^F{p) 



6 



)-— +6((p-—t + 7i 
M J V M 



(24) 



For another choice, say F2k{p) = (— l)^-F(p) and F2k+i{p) = 0, one obtains 



p(°)(0,p,t) = 7rF(p) 



6 



p 



5U 



p 



M 



t- 



71 

'2Ji 



(25) 



tH — 
M 2, 

All these are shown to be neutrally stable in Sec. El 

Equation ()23|) further indicates that there can exist higher-order multi-cluster 
phases (for C = 2,3,...) as well, if appropriate choices for Fk{p) are made. Recall 
again that the multi-cluster phase does not occur for stationary solutions since AW = 
for time-independent solutions. In other words, the multi-cluster must rotate with the 
frequency higher as the number of clusters grows; this suggests that the multi-cluster 
with large i should be difficult to observe. 



4. Stability of Stationary States 

In the previous work we have already shown that the stability of the incoherent 
phase depends on the solutions of the FPE, providing a plausible explanation as to the 
physical origin of the quasi-stationarity. We now extend the analysis further to include 
the case of the coherent phase. For this purpose, we write the FPE, setting A^^^ = A 
and 6 1 = 0, in the form 



dP V dP dP d 

dt M dcj) 'dp dp 



M " 



dp 



To probe the stability, we add a small perturbation to write 

P{(P,p,t) = Poi(P,p,t) + f{^,p,t) 
and accordingly 

A(t)=Ao(t) + Ai(t) 

dp I ^e'^^~'^ [Po(0,P,t) + /(0,P,t)]. 
Substituting these into ^lE^ . one obtains, to the lowest order. 



(26) 



(27) 



(28) 



dt 



p df 



M 

+ JAq sin 



+ JAi sin(0 — 



^dPo 
dp 



dp dp \M dp) ' 



(29) 



Since f{(f),p, t) and Ai(t) are periodic in 0, one can Fourier expand them in plane waves: 
f{<P,p,t) = J2[ due'^"^--'^f\{p,u) (30) 



Stability in a system of globally coupled rotors 



9 



and 



due-'^' [ dpf_,{p,uj), (31) 



where the integration over has been performed. Note here that the perturbed order 
parameter is proportional only to /_i(p, u;) (or to f+i{p,uj) if the order parameter 
has been defined to be A = {e~'''^'^~^^)) . Inserting these expressions into Eq. (j^^ 
and collecting coefficients of e^i^^l'-'^'t) ^ one finds the relations satisfied by the Fourier 
coefficients fk{p,u!). 

In the case of ferromagnetic coupling, the coherent phase (Aq 7^ 0) arises at 
temperatures below Tc, regardless of the presence of damping. The stationary solution 
in Eq. (fT^ can be written, with the help of Eqs. (fTB|) . (fTTj) . and ^T^ . in the form 



00 J 

Po{<P,p) = fM{p) f 



^ Jn(<t>^e) 



-00 



= fM{p) J2 e*"^'^"'^ (32) 

n=— 00 

where /m(p) = {2nMT)-'^/^ exp{-p^/2MT) is the Maxwell distribution and x = JAq/T 
as before. When x = 0, the above equation simply reduces to the Maxwell distribution, 
which is stable at temperatures above Tc. Our concern now is how the coherent phase 
gets its stability as the temperature is lowered below the critical temperature. Putting 
Eqs. and (jH^ into Eq. (j^ . one obtains the following equation for the Fourier 

coefficients fkip.oj): 

= ^[A^'-'\x) - A(^+i)(x)]/l,(p) J dp'l,. (33) 

We note here that the emergence of coherence contributes to the off-diagonal term 
in Eq. ()33|) and to the appearance of higher-order generalized order parameters, making 
the stability analysis non-trivial. When uj — kp/M = 0, we have a continuous spectrum, 
and for r = 0, Eq. (jH^I) becomes 



JAo d ~ J 



fVi = ^A^'~'\x)fM 



dp'l 



(34) 



2 dp 2 

It is easy to show by direct substitution that this equation has a solution of the form 

A(P.-) = I (35) 
I U, otherwise. 

It is of interest to note this is also the solution for F 7^ as well, since the term including 
F vanishes for the Maxwell distribution. For uo — kp/M 7^ 0, we have a discrete spectrum 
and may not solve the equation for the general case. Still we may proceed further if 
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we take the phase-only perturbation, namely, f{p,(f>,t) = fM{p)h{4),t). We then have 
the Fourier coefficient fk{p,uj) = fM{p)hk{oj) with hk{uj) being the Fourier coefficient 
of h{(j),t), which in turn gives {d/dp)fk{p,uj) = f'j^{p)hk{uj) and / dpfk{p,uj) = hk{uj). 
Further, the term including F vanishes identically in this case. Dividing Eq. (jH^ by 
uj — kp/M and integrating over p, we obtain 

/ifcM = - [A('^-i)(x) - A(^'+i)(a;)]x.H/^-i(a7) 

+ 2AoXfc(t^)[/ifc+i(t^) - hk-i{uj)], (36) 
where we have introduced the A;-dependent response function 



i I dp- 



'oj + kp/M ^'^'^^ 
and used Eq. ()A4|) . Some properties of this response function, which is frequency- 
dependent, are discussed separately in the Appendix. For a; 7^ 0, the recursion relation 
for the modified Bessel functions [TH] : 

(38) 



2k 

h+i{x) = — Ik{x) 

X 



Ik-i{x) 

leads Eq. to take the form 

2k 

hk - 2AoXk{uj){hk+i - hk-i) = —A^''\x)xkiuj)h^i, (39) 

X 

which needs to be solved. For k = 0, from Eqs. (j36|) and ()A2|) . we find Hq = 0, implying 

the absence of a constant term in the perturbation. Noting A^'^^x) = A^^'^^x) and 
Xk{i^) = —X-k{^^), we write the difference equation in the matrix form: 



A 



with the matrix 



A 



h.2 

V ••• / 

/ 1 







(40) 



2/x)AWxi 
A0X2 + (4/a;)A(2)x2 
(6/x)A(3);^3 



-AoXi 
1 

A0X3 





-A0X2 
1 



V 



(41) 



Here we have included the terms with only negative k values, reflecting that all 
order parameters are defined by Eq. ()18|1 . In order to have non-trivial solutions for 
h = (/i-i, /i-2, h-s, ...), one should have the vanishing determinant: 

e{tu) = detA = 0. (42) 

Let us ffist consider the limit Aq ^ or x = JAq/T —>■ 0, corresponding to the 
incoherent phase. In this limit all the off-diagonal terms vanish, since Io{x) 1 and 



(x) ^ (x/2)" so that A(") 

( l + Xi 
1 
1 

V 



(x/2)". Equation (pUj) then becomes 
\ ( h-i\ 



h-2 

V ••• / 



0, 



(43) 
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which leads to 



1 + xiiuj) = 1 + xiuj) = 1 + = 



(44) 



for non-vanishing h_i, while all other /I's are zero. The detailed analytic properties 
of the reduced response function X{^) = (2/^M)x(tt;), with the complex frequency 
u = Ur + iuJi, are presented in the Appendix. Equation (ji^ describes the condition 
for the incoherent phase with the Maxwell distribution, which is stable/unstable 
above/below ^H]- For a; 7^ 0, in principle we have to solve Eq. (j^^ including 
all the terms in Eq. (j4H) . Since this is very formidable, we instead consider just a 
few terms to explore how the stability of the solution changes. To this end, we write 
the determinant obtained when the first m Fourier components are 
kept. With only the first Fourier component h-i considered, Eq. (jl^ obtains the form 



uj) = l + -AWxi(^) 

X 

2T 

= 1 + — X(^) 

= 1 + TMx{uj) = 0. (45) 

for which Eqs. ()A11|) to ()A13|) yield cjj = as the only solution. Comparison of Eq. 
(I45|l with Eq. (jl^ shows the correspondence T = J/2 = Tc] this indicates that the 
solution is neutrally stable at the critical temperature, below which coherence develops. 
Including the next component h_2, one has 



Ao(Ao + -A(2))xi(c^)X2(c^) 

X 



0, 



which, with A*^-"^^ 



Ao and Xk{(^ 
1 



rp2 



x{!jj/k)/k, becomes 



X + Ax 



x{u:)x{uj/2) = 0. 



(46) 



(47) 



Since xi}^) and x(ct;/2) have the same pole structure, the real and the imaginary parts 
of Eq. (jUj) read 



Ree^'^\uo) = 1 + TRex(c^) + 



y2 



9 - h{x) 
x" + 4x-^ 



Ime^^^i 



X [Rex(^)Rex(^/2) - Imx(u;)Imx(^/2)] 



rp2 

= Tlmxi^j) H 



X + Ax 



h(x) 



0. 



(4^ 



(49) 



X [lmx{uj)Rex{uj/2) + Rex{uj)lmx{uj/2)] 

In the Appendix it is shown that tu^ = is a solution of Imx(u;) = 0, implying that this 
is also a solution of lme^'^\uj) = 0. For cuj > 0, Eq. (jl8|l becomes 

1 " 



fiy) - 1 



9 hix) 
x^ + 4x-^ 



/(y)/(2//2) 



(50) 



with y = uji^y M/2T. As y increases from zero to arbitrarily large values, the left-hand 
side of Eq. (|50|) decreases monotonically from zero to —1 while the right-hand side is 
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positive-definite for y > 0. This suggests tliat tliere is no solution for cjj > to make 
the system unstable. For tUj = 0, which corresponds to the neutral stability, we have 
Rex(i^) = — 1/T and thus Eq. (PH|1 reads 



X + 4x- 



0, 



hix) 

which leads to x = for the critical case. For Ui 
becomes stable, Eq. PHj) takes the form 



9iy) - 1 = o 



9 hix) 



hix) 



9iy)9iy/'2). 



(51) 

— \uJi\ < 0, for which the system 

(52) 



As shown in the Appendix, g{y) is a monotonically increasing function of y from unity 
to infinity, and accordingly, g(y)g{y/2) is also a monotonically increasing function of y 
in the same domain. Since the left-hand side of the above equation is monotonically 
increasing from zero to arbitrarily large values, Eq. allows a solution only for some 
range of x values. We have determined numerically the range of x values, in which there 
exits a solution for ?/ > 0, to find 



X 



T 



< X 



(2) 



1.32. 



(53) 



This indicates that the coherent solution is stable only at temperatures above Tq, at 
which Ao(T) and XcT/J meet. We see that the stable region does not extend to the 
zero temperature, presumably because we have included only the second component in 
our analysis (the first component is trivial). This may be resolved if one include higher 
components. Adding the third component hs leads to the following equation 



e^^Hio) 







24 ^ 



1 + Txioo) + 3T:^x(^ 
hix) 



x(^/2)x(^/3) 

(54) 



from which one can perform the similar analysis to find 

1] 



X 



1 + -f(y/2)fiy/3) 



2 , A hix) 
x^ + Ax ^ ^ 



for ujj > and 



_ 1 

~ 8 

1 + ^g(y/2)giy/3) 
hix) 



Ifiv) - 

hix) 
hix) hix) 



fivm 



fiy)fiym 



(55) 



X 



4x 



[9iy) - 1] 

hix) 



hix) hix) 



9iym 



9iy)9iym 



(56) 



for Ui < 0. Again, Eq. ()55|) does not have a solution for positive y since the left-hand side 
is less than zero while the right-hand side is greater than zero. Equation ()56|) is found 

(3) 

to have a solution for x < Xc ~ 1.51. Note here that Xc is increased substantially once 
the third component is included, implying that Tq, above which the coherent solution is 
stable, is lowered. We have performed this analysis, including up to four components. 
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and confirmed that this trend persists; this suggests the plausible conjecture that the 
coherent solution is stable down to zero temperature if all the Fourier components are 
included. 

We now turn our attention to the stability of the antiferromagnetic system for 
which there is no equilibrium order (Aq = or x = 0). The stability equation reads, for 
J=-\J\, 

1 - ^X(^) = 0, (57) 
which, depending on the sign of Ui (with Ur = 0), becomes 
l + (|J|/2T)/(y) = for.;, >0 

1 + |J|/2T = fortUi = (58) 

l + (|J|/2r)^(y) = ioTUJi<0. 

None of these equations has a solution, since f{y) and g{y) is positive-definite. This 
means that the antiferromagnetic system cannot have self-sustained deviation in the 
absence of the perturbation with a discrete spectrum. On the other hand, with the 
continuous spectrum u = Ur = kp/M, the system is neutrally stable at all temperatures. 



5. Stability of Non- Stationary States 



As mentioned in Section IHl the non-stationary solution exists only in the microcanonical 
ensemble (F = 0) with a continuous frequency spectrum, which can develop 
spontaneously. Our concern in this section is the stability of this solution, especially 
in the case of the antiferromagnetic interaction. Equation for stability reads, with 
Ao = F = 0, 



df P df . dPo 

—— = — — h JAi sm0— - — . 

dt Md<p ^ dp' 



(59) 



where Pq is given by Eq. (|22p. The last term in the above equation obtains the form 

dPo J 



JAi sin I 



dp 



- I dujc-''^' 
li 



dp'f-i{p',uj) 



X 



^-^ op 



exp ( -^^M ^kip) 



J_ 

2i 



$:||cio.[(e^('=-^)^-e^(-^)^)]exp 



kp. 



xFu{p) / dp'U{p\u) 



J 

2i 



/ due 



i{k(j>~LL)t) 



d_ 
dp 



Fk-i{p,Lj) - Fk+i{p,uj) 



with 



F,{p,u;)^F,{p) I dp'l,{p',uj-^] 



(60) 



(61) 
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which leads to the equation for the Fourier coefficients: 



UJ 



kp 
M 



J_d_ 
2dp 



Fk-.i{p,uj) - Fk+i{p,uj) 



(62) 



Since we are deahng with the perturbation of the non-stationary state with a continuous 
spectrum, the frequency of the perturbation should satisfy a; — kp/M ^ 0; otherwise, 
there would be no perturbation at all. This allows us to divide Eq. hy uj — kp/M 
and to integrate over p. For A; = — 1, we have 



dp'f-i{p',uj) = ^ 
d 



dp 



X 



F_2(p,u;) - Fo{p,Lj) 



(63) 



while for k ^ —1, f^ij)-, t^) is determined by f-i{p', u!±kp/M) through Eqs. (jF)T|l and (jHS))- 
It is thus enough to have non- vanishing /_i(p', u). Now suppose that u = UqIS a solution 
of Eq. i.e., f dp f^i{p',uj) 7^ for = ujq. If we write J c/p'/_i(p', cu) = K5{ijJ — ujq), 

then J dp'f-i{p', on + ^) = K6{uj + |f — t^o)- Integration over u gives 



JM 



dp 



nip) 

p + MiOo 



JM 
JM 



dp 
dp 



F-2{p) 
{p - Muo) 

p — Mujq ' 



(64) 



where the last line is obtained by integration by parts. Hence the frequency of a self- 
sustained oscillation and accordingly the stability is, similarly to the stationary case 
[Eq. (HH)], determined by 



1 + 



JM 



dp 



nip) 



F'-2{P) 



0. 



(65) 



p + Muq p — Mljq 
The stability condition is thus entirely the same as that of the stationary case 
except that we now have two momentum distributions: From Eqs. ()A8|) and ()A9|) with 

Mujq = ujr + oJi, we have 



+ 



JM 

IZodP 



dp 



{p + COr)Fl,{p) {p-Ur)FLM 



nip) 



{p + Ur)"^ + Luf 



iP 



-.2 



+ UJ. 







(66) 



0, 



{p + UJrY + Uj1 {p - UJ^)"^ + Uj"^ 

for uoi > 0, for which the system is unstable as the perturbation grows in time. In 
the opposite case {ui < 0), the perturbation dies out to make the system stable. The 
condition for this is given by 



JM 



+ 27r 

00 

dp 

00 

+ 27r 



dp 



ImF^ 



{p + CO r)Fl,{p) {p-Ur)F'_M 

{p + u, 
i-u) 

F'oip) 



Co) - ImF' . 



{p - UJr) + Cj 

-Cj)] = 



FUip) 



_ip + UJr 



ReF^i 



-UJ 



{p-UJr)^ + ul 

ReFL2i-Co)] = 0. 



(67) 
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Finally, in the neutral case (cjj = 0), the condition simply reads 



^ V I dp 



JM 



(68) 



P + UJr p — UJr 
Fl,{-Ur)-F'_^{-Ur)=Q, 

where V stands for the principal part. Our next task is to determine stability for 
specific distributions of Fq{p) and F_2(p). Since most dynamical calculations, for 
both ferromagnetic and antiferromagnetic system, have used the so-called water-bag 
distribution, we also consider the momenta to be distributed uniformly in the range 
[—a, a]: 

Fo(p) = ±F_2(p) = ^ (69) 

Substitution of Fq{p) = ±F^2{p) = i2a)-^[S{p+a) - S{p-a)] into Eqs. 
depending on the sign of Ui, determines the frequency Uq. 

We first consider the case Fo{p) = F_2(p). From the second equations of Eqs. (jHEl) 
(for Ui > 0) and (jU7j) (for cUj < 0), we find Ur = 0, while there is no solution to satisfy 
the first ones. When uJi = 0, again there is no solution to satisfy the first equation 
of Eq. ()(j8|l . This indicates that there is no self- sustained oscillation in the system. 
Note, however, that the system is neutrally stable as it has a continuous spectrum 
{u = kp/M). Next, when Fo{p) = — F_2(p), one finds 



" ^^/j4~ ijd) ' = ° for a < Or (70) 



uJr= ±\l + , = ioT a>aR (71) 



with = V JM. In the microcanonical ensemble one may relate the average 
kinetic energy with the temperature: T/2 = (jP')/2M = a'^/QM, from which one has 
Tr = aj^/3M = J/3 [in]. Accordingly, it is concluded in this case that the rotating 
solution is neutrally stable for T > Tr and becomes unstable below Tr. Note here that 
Tr is lower than the equilibrium critical temperature Tc = J/2. 

For the antiferromagnetic system (J < 0), we replace J = —\J\ in Eq. (fTOj) to 
obtain, for ujj = 0, 



while there is no solution for uJr = 0- We therefore conclude that the antiferromagnetic 
system is neutrally stable for all a, i.e., at all temperatures. In Sec. 01 we have shown 
that the bi-cluster state is allowed by the rotating distribution. The result obtained here 
that this non-stationary solution is neutrally stable at all temperatures thus suggests 
an alternative explanation as to the origin of the spontaneously formed bi-cluster state 
in numerical simulations, which retains its form for quite a long time ^21- This keeps 
parallel with the emergence of quasi-stationarity in the ferromagnetic system, associated 
with the neutral stability [T!^ . 
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6. Conclusion 

In this paper, we have presented a detailed analysis of the system of globally coupled 
rotors. Starting from a set of Langevin equations and their corresponding FPE, which 
includes the micro canonical ensemble approach as a limiting case, we have found a class 
of solutions and studied their stability. The standard canonical distribution constitutes 
a simultaneous solution of the canonical and the microcanonical ensembles, and thus 
describes the same equilibrium behavior in both ensembles, leading to the coherent 
phase (characterized by a nonzero mono-cluster order parameter, i.e., A^^^ 7^ 0) below 
the critical temperature Tc in the ferromagnetic system. The stability of the coherent 
phase is governed by an infinite-order difference equation, the behavior of which may 
be understood by considering successively higher-order terms (i.e., Fourier components 
in the perturbation). It has been found that the coherent phase is stable above some 
temperature Tq, which is finite if one includes only a few lowest Fourier components. 
As more components are considered, Tq tends to decrease toward zero; this leads us 
to surmise that the infinite number of Fourier components would stabilize the coherent 
phase down to zero temperature. Namely, it is expected that the stability equation, if 
treated exactly, leads to the stability of the coherent phase at all temperatures below 
T 

c- 

We find the more interesting possibility for the non-stationary (rotating) solution 
with regard to dynamical order. Dynamical order, manifested by multi-cluster 
motion, is allowed for both ferromagnetic and antiferromagnetic interactions. Unlike 
a ferromagnetic system, in which dynamical order ceases to exist below the temperature 
Tr, dynamical order is observed to be neutrally stable down to zero temperature 
in the antiferromagnetic system. This suggests an alternative explanation as to 
the origin of the spontaneous formation of the bi-cluster phase in the system of 
antiferromagnetically coupled rotors. This is in parallel with the explanation that the 
quasi- st at ionarity observed in ferromagnetically coupled rotors is related to the neutral 
stability of the stationary solution in the incoherent phase below the equilibrium critical 
temperature [T3j. 

To conclude, we have introduced a unified approach for both the canonical ensemble 
and the microcanonical ensemble, based on the Fokker-Planck equation. Depending on 
the ensemble, the Fokker-Planck equation admits a few solutions which have implications 
on some remarkable features (quasi-stationarity in ferromagnetic systems and bi-cluster 
motion in antiferromagnetic systems) observed in numerical experiments. We provide 
natural explanations for the origin of those seemingly unrelated features within the same 
context. Finally, we point out that our approach is based on an effectively one-particle 
dynamics, exact for an infinite number of particles and does not reflect the instabilities 
that may be caused by the finiteness of the number of particles. 
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Appendix A. Properties of Xk{^^) 

In this appendix we describe some properties of the response function Xki^^) for the 
Maxwell distribution fuip)'- 

First, for /c = 0, we have 

Xo(-) = ^ / dpf-M^ = (A2) 

since fM^p) is an odd function. We next write k —k and change the integration 
variable p to —p in Eq. ()A1|) to get 



X^k{uJ) = - dp- 



2 J oj- kp/M 



2 J ' "^'u + kp/M 

again noting that f'^ip) is an odd function. Similarly, it is straightforward to show that 
Xki-^) = -X-ki^) = Xki^^)- (A4) 
Further, Eq. ()A1|) can also be written as 



Xfc(^) = / dp- 



where 



2k J p+Mu/k 
= Ixi^/k), (A5) 

is the response function already defined in Sec. |3 Although we consider here the 
Maxwell distribution, the properties given above hold for any momentum distribution 
/o(p), only if it is an even function of p. We now proceed to evaluate this function. 
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paying attention to the simple pole at p = —Mu on the complex p-plane. Setting 
Muj = u and making analytic continuation uj = uJr + iuJi, we obtain x{^) i^i the form 



JM 



-UJ) 

-Co) 



With the tilde sign omitted for convenience, the real part reads 



Rexiuj) 



for cijj > 
for Ui = 
for Ui < 0. 

for tUj > 

for Ui = 
for cu,- < 



(A7) 



while the imaginary part is given by 



Imx(c^) 



IZo dp 



flfip) 



for > 
for uji = Q (A9) 
^ + 27rRe/j^j(— Lu) for tUj < 0. 



We next write 



{p + iOr) + iO.^ 



+2 I LOi COS 



UJr COS , ,„ + Ui sm ■ 



MT 



UJr sm ■ 



MT 

UJfUJi \ 



(AlO) 



MT "'~"'MTJ. 
= Reflj{uJr+iuJi) +ilmflj{uJr+iuJi), 

from which it is obvious that Kef'j^{ujr + iuJi) = for cj^ = and lmf'j^j{iuji) = 
-(27rM3T^)-i/2^.e'^?/2MT_ ^Ve thus conclude that u;^ = is a solution of Imx(c^) = 0, 
since fj^ip) is an odd function of p, which makes the integrals vanish in Eq. ()A9|) . We 
now evaluate the integral of Rex(co'). For uJr > 0, the first equation in Eq. ()A8|) becomes 

m 

1 



Rexiuj) = - -[1 - Vnye^ erfc(?/)] 



T" 



(All) 



''"dt 



(A12) 



with the scaled variable y = uji^ M j^T ^ where 
erfc(?/) = ^ 

is the complimentary error function. For uji = ~\uJi\ < 0, it is straightforward to show 
that the last equation in Eq. ()A8|) becomes 

1 



Rexiuj) 



y\ey (2-erfc(|y| 



Stability in a system of globally coupled rotors 



19 



^-f9iy)- (A13) 

For LJi = 0, we have Imfj^j^oUr + ioJi) = and the second equation in Eq. ()A8|) simply 
reduces to Rex(i^) = — 1/T. Note that f{y) is a monotonically decreasing function of 
y, varying from unity to zero as y grows from zero to arbitrarily large values. On the 
other hand, g{y) increases monotonically with y, from unity to arbitrarily large values. 
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